Unconventional superconducting states induced in a ferromagnet by a rf-wave superconductor 
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We develop a quasi-classical theory for the superconducting proximity effect in a ballistic ferromagnetic layer 
in contact with a d-wave superconductor. In agreement with recent experiments we find that the density of states 
oscillates around the normal state value with varying the thickness of the ferromagnetic layer We show that the 
phase, the amplitude, and the period of these oscillations depend on the orientation of the superconductor This 
effect reveals spatial oscillations and anisotropy of the induced superconducting correlations in the ferromagnet. 



Proximity effects in hybrid structures of superconductors 
and ferromagnets provide the possibility for the controlled 
studies of the coexistence of ferromagnetism and supercon- 
ductivity. Superconducting correlations penetrate into normal 
metals by means of a special scattering process at the normal 
metal-superconductor (NS) interface known as Andreev re- 
flection yJ]: an electron propagating in the normal metal can 
cross the NS boundary by reflecting as a hole and transfer- 
ring a Copper pair into the superconductor. The Andreev re- 
flected hole is correlated with the incident electron and carries 
information about the phase of the order parameter of the su- 
perconductor. The existence of correlated electron-hole pairs 
corresponds to a non-vanishing superconducting pair ampli- 
tude (order parameter), which gives to the normal metal the 
typical characteristics of the superconducting states |2]. 

The peculiarity of this proximity effect in ferromagnetic 
metals comes from the fact that the Andreev reflection from a 
singlet pairing superconductor is accompanied by an inversion 
of the spin direction, and consequently the correlated electron- 
hole occupy opposite spin bands. The spin splitting exchange 
field of the ferromagnet causes a phase shift of the correlated 
Andreev electron-hole pairs, which results in a spatially os- 
cillatory behavior of the induced pair amplitude fj]. Such an 
inhomogeneous superconducting state coexisting with ferro- 
magnetism is strikingly similar to an FFLO state in a bulk su- 
perconductor with spin splitting originally predicted by Fulde 
and Ferrell and by Larkin and Ovchinnikov |4]. 

One of the most interesting manifestations of the oscil- 
lations in the proximity pair amplitude is the oscillatory 
behavior of the density of states (DOS) in ferromagnet- 
superconductor bilayers as function of the ferromagnetic layer 
thickness. Kontos et al |5] have observed this effect exper- 
imentally in thin ferromagnetic layers connected to a con- 
ventional s-wave superconductor. Theoretically, the super- 
conducting proximity effect in ferromagnetic layers has been 
studied extensively |3, 6]. It has been shown |7| that the the- 
ory based on the quasi-classical formalism |8] is in a quanti- 
tative agreement with the experiments. 

In parallel, there also has been much attention to the prox- 
imity of a ferromagnet to high-Tc superconductors, which are 
widely believed to have a dominant d-wave pairing symmetry 
10]. Most of the studies have concentrated on the influence of 
the exchange field on the zero bias conductance peak (ZBCP) 
as the most remarkable feature in the tunneling spectroscopy 



of a junction between a high-Tc superconductor and a normal 
metal ZBCP is the result of a zero-energy Andreev 

peak (ZEAP) in the DOS, which originates from change in 
the sign of the d-wave order parameter under a 7r/2 rotation. 
Many effects including a splitting of the ZBCP and the devel- 
oping of a zero bias conductance dip (ZBCD) have been found 
and attributed to the exchange splitt ing 1 1 2ll 1 3ll 1 411 . 

Most recently, Freamat and Ng \ly\ have performed tun- 
neling spectroscopy measurement on multilayered junctions 
of d-wave superconductor and ferromagnets. They observed 
an oscillatory behavior in the tunneling conductance spectra 
with the thickness of the ferromagnetic domain. This effect 
has been considered as the d-wave analog of the effect ob- 
served in |5] and explained as the signature for an FFLO state 
induced in the ferromagnetic layers by the d-wave supercon- 
ducting layer. However one would expect that the anisotropy 
and the sign change of the d-wave superconducting order pa- 
rameter leads to an anisotropy and phase effects in the os- 
cillation of the DOS, which distinguish it from the case of a 
conventional s-wave superconductor The goal of the present 
work is to study this effect theoretically. 

In this letter we study the proximity effect at the interface 
between a ferromagnetic metal and a d-wave superconductor 
We model the ferromagnet as a thin ballistic film with rough 
boundaries which is connected to the superconductor through 
a disordered interface of finite transparency. We make use 
of the quasi-classical formalism |8] to calculate the density 
of states in the ferromagnetic layer In correspondence with 
the recent experiments, we find that at sufficiently strong ex- 
change field the DOS oscillates as a function of the thickness 
of the ferromagnetic layer We show that the phase, the am- 
plitude and the period of the DOS oscillations are affected by 
changing the orientation of the d-wave order parameter This 
effect, which reveals the highly anisotropic nature of the in- 
duced superconducting correlations in the ferromagnet, can 
also be used as a further test of the d-wave scenario in high-Tc 
superconductors. 

We consider the proximity system shown in Fig. Q A thin 
ferromagnetic layer (F) of thickness d is contacted by a d- 
wave superconductor (S) on one side and covered on the other 
side by an insulator. We consider a clean structure at both 
of the S and F sides of the FS-interface. At the same time 
the FS-interface itself may contain disorder and band mis- 
match which enhance backscattering of quasi-particles from 
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with spin a (= ±1) is described by an independent equation 
for the corresponding matrix Green's function gcr(^', vp, r), 
which reads 



S 




FIG. 1: Sketch of a ferromagnetic layer (F) in connection with a 
d-wave superconductor (S) . A typical classical trajectory is also 
shown. It starts at S and extends into F making an angle Q\ with the 
interface normal and then, after several reflections from the insulator 
and FS-interface, returns into S at an angle ^2- 

this boundary. F is characterized by a mean exchange field h 
which is assumed to be constant in the F-layer and vanishes 
in S. S is characterized by a dr^i^yi-^NW^ order parameter of 
the form Ac((6') — A cos [2(0 — %)], where Q is the angle that 
the Fermi velocity Vp makes with the normal to the interface 
and X is the angle between the crystallographic a axis of S 
and the interface normal. For calculation we adopt the Eilen- 
berger equations |8|] in the clean limit for this system. In the 
absence of spin-flip scattering in F, transport of quasi-particles 



-ivp\/ga{E, Vf, r) = [{E + ah{r))f3 
-if2A(r,VF),5.(£;,VF,r)]. (1) 

Here denote the Pauli matrices and A(r,VF) is the super- 
conducting order parameter, which is taken to be real. We 
neglect the spatial variation of the order parameter close to 
the FS-interface and take A(r, vf) = Ad(vF) in all points of 
S. 

Equation Q can be conveniently solved along a classical 
trajectory. A typical trajectory is shown in Fig. ^ It starts 
at S and extends into the F layer making an angle di with 
the normal to the FS-interface. After several reflections from 
the insulator and the FS interface, which makes a path of to- 
tal length I in F, the trajectory returns into S at an angle 02- 
The superconducting order parameter that an electron feels in 
the beginning and at the end of the trajectory are respectively 
Ad{Oi) and Ad{92)- We have solved Eq. ([0 and found that 
the DOS on a given trajectory only depends on di, 6*2, and I 
and given by 



7V(£;, 01,02,0 = ^ Yl Re[trf3g.(£^,VF,r)] = ^ ^ Re{-i tan [M/«f + + «(02) + 77(^1, ^2)/2]}, (2) 
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where No is DOS at the Fermi level in the normal state, 
a{e) = arcsin[i;/|Ad(6')|]/2, = {E + ah)/vF and 
r]{0i,e2) = 7r[l-sign(A<j(0i)Ad(02))]/2 is the tt phase shift 
resulting from the possible change in the sign of the order pa- 
rameter at the beginning and the end of the trajectory. 

The total density of states N{E) is obtained by averaging 
the DOS per trajectory Eq. (|3 over all the possible values of 
01,02, and /: 

N{E)^ j dld0xd02p{l,0i,e2)N{E,l, 01,02), (3) 

where p{l, 0i,02) is the distribution function of the trajecto- 
ries. 

To determine p{l, 0i, 02), we model F as a weakly disor- 
dered thin layer bounded by a rough surface and a rough 
FS-interface with a constant transparency t. We introduce 
g{0i,02) as the coiTelation function between the incoming 
and outgoing directions 0i and 02- This function determines 
the specularity of the scattering from the boundaries. For 
perfectly flat boundaries with purely specular scattering, 0i 
and 02 are completely coiTelated and g S{0i +02)- For 
the boundaries with strong roughness the dominantly diffu- 
sive scattering destroys any correlation between 0i and 02, 



and g{0i, 02) is constant. In the general case of an arbitrary 
strength of the roughness we assume g = C^^exp(— (0i + 

02)'^/{Trz)^), where C = fj^^^ d0xd02g{0i,02) is a normal- 
ization factor and z measures the degree of the roughness. The 
limits of z <C 1 and z ^ \ correspond to weak and strong 
roughness respectively. As we will see below the parameter 
z play a crucial role to reveal the anisotropic nature of the 
induced superconducting correlations in the F-layer 

The full distribution function p[l, 0i, 02) can be written in 
terms of g{0i,02) as 



p{l 
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i,02)^g{0i,O2)mi 
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(4) 



where p{l) is the length distribution of the trajectories regard- 
less of the values of 0i and 02, and r — 1 — t. For t — 1 
the delta function in Eq. © assure that the total length of the 
trajectory is given hy I — d/ cos 0i + d/ cos 02- In the limit 
of small transparency of the FS-interface t ^ 1, p(0 can be 



FIG. 2: Effect of weak exchange fields on the density of states when 
the transparency of the FS-interface is small (here t — 0.1). Nor- 
malized DOS N{E)/No as a function of dth/v^ (dt ~ d/t) at (a) 
different energies when X = and z — 0.5, and (b) zero energy 
when X = and for different strength of the roughness of bound- 
aries z. The exchange field suppresses the typical superconducting 
features of the DOS (the zero energy Andreev peak when x = ''r/4, 
and also when x = if the boundaries are rough, and the vanishing 
DOS at _E = when the boundaries are specular and x = 0). 

approximated by the exponentially decaying function 0] 

p{l) = ie-'/'e(/ - 2d), (5) 

where I ~ 2d\n{£imp/d)/t — 2dt is the mean trajectory 
length in the limit of small d/iimp {^imp is the elastic mean 
free path). 

Combining Eqs. (|2}, Q, 0, and Q we find the following 
result for the total DOS 

N °° 
N{E)^^Y. I] [i + ?-p(2?^fc.d)]^, (6) 

o"— ±1 n— — oo 

where p{k) = e"^'^'"^ / {ikl + 1) is the Fourier transform of 
p{l) and 

Jn{k)= / deide2g{ei,e2)fniei)fn{02)e-'"''^'''''\ 

J-rr/2 

Jn{V) - < ^-i2nk„d/cos0-2\n\l3{e) f^j. |Ad(6')| ' 

Here 0(6*) = arccos [E /\/S.d{d)\]/2 is the Andreev phase and 
=acosh[i;/|Ad(0)|]/2. 
Equation (|6j expresses the total DOS in terms of known 
functions. In real samples the F-layer is expected to have a 
nonuniform thickness due to the large scale roughness of the 
boundaries. To take this into account we average (|6} over a 
Gaussian distribution of the thicknesses around a mean value 
d. This leads to a smearing of the sharp features in the DOS 
which arise from the lower cutoff in the length distribution. 
The qualitative behavior, however, will not be changed. In 
practice, we have taken the width of the distribution to be of 



-1 -0.5 0.5 1 
E/A 

FIG. 3: Oscillations of the proximity density of states with the F- 
layer thickness d. Dependence of the normalized DOS N{E)/No 
on d when the interface is highly transparent (t — 1) and h = 3A, 

X = 7r/4, z — 1. 



order 10%, which corresponds to the conditions of the exper- 
iments Bill. 

Let us start by analyzing the effect of a weak exchange field 
on the proximity DOS of a thin F-layer with d <^ vp/A and 
t ^ 1. For the orientation x = '^/'^ the main feature in 
the DOS of a normal layer (h — 0) is a sharp ZEAP for any 
strength of the roughness z. In Fig. the normalized DOS 
N{E)/Nq versus dth/vp is plotted for different energies and 
values X — i'/4, z = 0.5. The effect of the exchange field 
is to split the ZEAP into distinct peaks at finite energies. In- 
creasing h further leads to the shifting of the splitted peaks 
toward higher energies, which is associated with decreasing 
the height of these peaks. The DOS at _E = decreases with 
h and goes through a minimum as dth/v-p is increasing. Thus 
the exchange field can suppress the ZEAP in the DOS and 
develop a dip at ii^ = 0. This will lead to a ZBCD in the 
conductance spectra at low temperatures, which also has been 
observed in the experiments L1Z1. 

For the orientation x = the ZEAP is absent for perfectly 
specular boundaries (z = 0). The ZEAP, however, appears 
when z is finite. This is the result of diffusive scattering 
caused by the roughness at the boundaries, which allows for 
the sign change of the superconducting order parameter at two 
sides of a sizable fraction of trajectories even when x = 
d. In Fig. IS) we plot the normaHzed DOS Af(0)/iVo at 
E = Q versus dth/vp when x = and for different values of 
z. For the normal layer (h — 0) the height of the ZEAP in- 
creases with increasing z. The effect of the exchange field for 
a given value of z is to suppress the superconducting features 
(zero DOS for z = and ZEAP for finite z ) in the DOS. At 
higher exchange fields when dth/vf^l, the zero energy DOS 
for all values of z approach each other and becomes close to 
the normal state values A^o- 

At strong exchange fields when dh/vp > 1, the DOS at 
all energies oscillates around the normal state values iVo as a 
function of dh/v^. This oscillatory behavior of induced su- 
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FIG. 4: Effect of the d-wave symmetry on the oscillations of the 
DOS. The reduced DOS at zero energy [A'^(0)/A*'o — l]/t versus 
dh/vj: for different orientations x when the boundaries are (a) spec- 
ular (z — 0) and (b) diffusive (z = 1). In both cases the amplitude 
and the phase of the oscillations are modulated by varying x- For 
diffusive boundaries the period of the oscillations is also different for 
two cases of x ~ 7i"/4 and x ~ 0. (c) shows the transition from the 
specular to the diffusive case by varying z for fixed x = 7r/4 . 



X = 7r/4. We note that the period of the oscillations is almost 
the same for different x and equals to 7r/2, the period for the 
case where S has a s-wave order parameter 0]. 

The effect of an anisotropic order parameter is more pro- 
nounced for the case of diffusive boundaries. This is shown 
in Fig. 1^ where the DOS oscillations at = are pre- 
sented for the case of z = 1. Now, in addition to the variation 
of the phase and the amplitude, by changing the S orienta- 
tion from X — to X = 7r/4 the period of the oscillations 
is also changed. Fig. |^ shows how the change in the pe- 
riod occurs when the strength of the roughness z varies for 
X — 7r/4. We can understand these effects by noting that 
the spatially oscillating order parameter in the F-layer has a 
direction-dependent amplitude and phase resulting from the 
anisotropy and the sign change of the d-wave order parameter 
of S. The change in the phase, amplitude and the period of 
the DOS oscillations is the result of unconventional induced 
superconducting correlations 

In conclusion, we have investigated theoretically the su- 
perconducting proximity effect in a thin ferromagnetic layer 
in contact with a c?-wave superconductor In correspondence 
with recent experiments 1 15] we have found that at sufficiently 
strong exchange fields the density of states in the ferromagnet 
oscillates around the normal state value as a function of its 
thickness. The phase, the amplitude, and the period of the 
oscillations depend on the orientation of the superconductor 
This direction-dependence is the signature of an unconven- 
tional oscillatory superconducting state induced in the ferro- 
magnet by the proximity to the c?-wave superconductor. 

We thank M. Freamat, P. Fulde, Yu. V. Nazarov, and H. 
Schomerus for useful discussions. 



perconducting correlations has been observed in the experi- 
ments 1 15]. Fig.|2|represents these oscillations resulting from 
our calculations. The energy dependence of the normalized 
DOS N{E)/Nq is plotted for different values of the F-layer 
thickness when the FS-interface is highly transparent {t = 1), 
h — 3 A and x — tt/4:. We have taken z = 1, modeling rough 
boundaries. In this case the period of the DOS oscillations is 
roughly tt/A. 

Let us now analyze the effect of the c?-wave symmetry of 
the superconducting order parameter on the oscillations of the 
DOS in F. For this we examine the oscillatory behavior of the 
DOS at the Fermi level (E — 0) for different orientations of S. 
We consider the case of i ^ 1 where for h> vp/d the devia- 
tions of the DOS from the normal state value A^o are of order t. 
In Fig. Hi we have plotted the reduced DOS [iV(0) /Nq - l]/t 
versus dh/vp for the case of specular boundaries (z = 0) at 
different x- The variation of x affects the phase and the am- 
plitude of the oscillations. Varying x from to 7r/4, the phase 
is shifted by tt and the amplitude is decreased. The tt phase 
shift originates from the sign change of the order parameter, 
which occurs for all the trajectories in specular boundaries if 
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